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Abstract: This article brings in two new discrete distributions: multivariate Binomial distri- 
bution and multivariate Poisson distribution. Those distributions were created in eventology 
as more correct generalizations of Binomial and Poisson distributions. Accordingly to even- 
tology new laws take into account full distribution of events. Also, some properties and 
characteristics of these new multivariate discrete distributions are described. 
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1 Introduction 



Distribution of probabilities is one of principal idea in theory of probabilities and math- 
ematical statistics. Its determination is tantamount to definition of all related stochastic 
events. But trials' results extremely rarely are expressed by one number and more fre- 
quently by system of numbers, vector or function. It is said about multivariate distribution 
if some regularity is described by several stochastic quantities, which are specified on the 
same probabilistic space. Thereby, it is involving for representation of behavior of the 
random vector, which serves a description of stochastic events, more or less near to re- 
ality. This work came in connection with appeared scientific necessity of assignment two 
new multivariate discrete distributions, which are naturally following from eventological 
principles and taking into account specific notion from probability theory namely arbitrary 
dependence of events, random quantities, tests. 

2 Binomial multivariate distribution 

Let there are finite sequence of n independent stochastic experiments. In the result of 
experiment i can ensue or not events from N-set X^ of events £ X^. Eventological 
distributions of sets of events 3£( l \i = 1, ...,n agree with the same eventological distribu- 
tion {p(X),X C X} of certain iV-set X of events x £ X, which aren't changing between 
experiments. 

Such scheme of testing is called multivariate (eventological) scheme of Bernoulli testing 
with producing set of events X, and each of random quantities 

n 

i=i 

obey the Binomial distribution with parameters n, p x = P(x), while random vectoJl 
£ = (£x,x £ 3C) obey the Binomial multivariate (N-variate) distribution with parame- 
ters (n, {p(X), 
^ X C X}). 

Probabilities of Binomial multivariate distribution, which is generated by N-set of events 
X, are defined for every integer-valued collection h = (n x ,x £ X) £ [0,72,]^ by the formula 

fe ft (n;p(l),0^ ia)=P(( = n) = Pfe = n x ,x £ X) = ^m^n; {p(X),X C X}), 

h 

where 

mfl (n; {p(I),ia})=P(( = n)=p((((4ia) = (n(X),XCX)) = 

- n ^ <x) 

11 n(X)\ XQX 

xcx 

3 In eventology, in general, and at this context, in particular, the notion jj vector is using in extended 
sense as disordered finite set or disordered finite collection of some elements. 
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are probabilities of 2 Ar -variate Polynomial distribution of a random vector £ = (£(X), X C 
X) with parameters (n;{p(X),X C X}), which is generated by 2 N -set of terrace-events 
{ter(X), 

ICI}, which has biunique correspondence to the given Binomial multivariate distribu- 
tion; summation is made by all 2 7V -variate sets h = (n(X),X C I) e S 2 from 2 N - vertex 
simplex S 2 " , i.e. such as 

n= ^n(X), 

but which are meet the N equations 



^n(I),iG X. 



= 

2.1 Binomial one-variate distribution 

When N = 1 (i.e. generating set X = {x} is a monoplet of events) Binomial one-variate 
distribution of a random quantity £ x coincides with the classical Binomial distribution with 
parameters (n;p x ). In other words, probabilities of the Binomial one-variate distribution 
have classical format 

b nx (n;p x ) = P(£, = n x ) = C>^(1 - p x ) n - n *, 0^n x ^n. 

2.2 Binomial two-variate distribution 

When N = 2 (i.e. generating set X = {x, y} is duplet of events) Binomial two-variate 
distribution of a random vector £ = ((, x ,(,y) = (£x,x £ X) is defined by four parameters 
(n;p(x), 

p(y),p(xy)\ wherej 

p(x) = P(x fl y c ), p{y) = P(x c fl y), p(xy) = P(x fl y). 

Probabilities of Binomial two-variate distribution are calculating for any integer- valued 
vector n = (n x ,n y ) G [0,n] 2 by the formula 

bn(n;p(x),p(y),p(xy)) = P(| = h) = P(£ x = n x , f„ = n y ) = 

min{n x ,n y } 

= ^2 r nn(n;p(0),p(x),p(y),p(xy)), 

n(xy)=ma,x{0,n x +n y —n} 



where 



mn(n;p(0),p(x),p{y),p{xy)) = P(| = n) = 
P((Z(0),Z(x),Z(y),Z(xy)) = (n(0),n(x),n(y),n(xy)) 

- [ P (0)] n{0) ip(x)} n(x) [p(y)] n(v) Wv)\ <xv) 



n ! 



n(0)\n(x)\n(y)\n(xy)\ 



Obviously, p(0) = 1 — p{x) — p{y) — p{xy). Hereinafter are used next denominations: p x = ~P{x) 

.2 

y 



p(x) + p{xy), p y = P(y) =p(y) +p(xy), Kov xy = p(xy) - p x p y , o\ = p x (l -p x ), o y = p y {l - p y ) 



are probabilities of the Polynomial 4-variate distribution of random vector £ = (£(0), £(x), 
£(y) , £(xy)) with parameters (n; p{0) , p(x) , p(y) , p{xy))\ summation is made by all sets 
ft = (n(0) , n(x) , n(y) , n{xyj) such as n = n(0) + n(x) + n(y) + n(xy), for which are true 
2 equations as well n x = n(x) + n(x,y), n y = n(y) + n(x,y), and can be turned into 
summation by the once parameter n(x, y) within Frechet bounds, since when n x and n y 
are fixed, then all quantities n(0), n(x), n(y), n(xy) can be expressed by one parameter, for 
instance, n(xy): 

n(x) = n x — n(xy), n(y) = n y — n(xy), n{0) = n — n x — n y + n(xy), 
which is varying within Frechet bounds: 

max{0, n x + n y — n} ^ n(xy) ^ min{n x , n y }. 
Also, formula can be written in the next manner: 

bn(n;p(x),p(y),p(xy)) = P(| = n) = 

min{n x ,n y } 

= b(0)r [r(x)r [r(y)P E ' G n n ^\h) [t(x, y)f xy) , 

n(xy)=max{0,n x +n y — n} 

where 

n\ 



in — n x — n y + n{xy))\{n x — n{xy))\{n y — n(xy))\n(xy)\ 
is two-variate Binomial coefficient, and 

p{x) p(y) p(0)p(xy) 

t(x) = — — r, T(y) = — — -, r(x,y) = — ; - - - 
1 J p(0) K } P(0) 1 ^ P(x)p(y) 

is first and second degree multicov aviations of events x and y. 

Vector of mathematical expectations for Binomial two-variate random vector (£ x ,£y) 
is equal to (E£ X ,E£ X ) = (np x ,np y ), and its covariation matrix can be expressed through 
covariation matrix of random vector (l x , l y ) of indicators of events from the generating set 
X = {x, y} and is equal to 

np x (l-p x ) nKov xy \_ n fpx(l-Px) Kov X2/ 



nKov xy npy(l-p y )J \ Kov xy p y (l-p y ) / 

Covariation matrix of the centered and normalized Binomial two-variate random vector 

'tx-np x iy-npy 



a x Oy 



is expressed through covariation matrix of the random vector 

Px Py 



0~x 0~y 



centered and normalized indicators of events from X = {x, y} and is equal to 

n np xy \ _ n ( 1 p xy 



/TLPxy 71 J \Pxy 1 

tion coe 

indicators of events from X = {x,y}). 



where p xy = °° Y * y is correlation coefficient for random quantities l x and l y (i.e. for 
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2.3 Characteristics of Binomial multivariate distribution 

Vector of mathematical expectations for multivariate random vector (£ x ,x G X) is equal to 
(E£ x , x G X) = (np x , x G X), its covariation matrix is expressed through covariation matrix 
of random vector (l x , x G X) of indicators of events from the generating set X and is equal 

to 

na 2 x ... nKov^X / a* ... Kov xy 

K nKov xy . . . na^ J yKov^ ... a 2 

where a 2 x = p x (l - p x ), Kov xy = -p x p y when x^y. 

Covariation matrix of generated by partition centered and normalized Binomial multi- 
variate random vector 

'£x - np x ( c ^ 

is expressed through covariation matrix of random vector 

(l x _-p x )_^ x G X 

0~x 

of centered and normalized indicators of events from X and is equal to 

' no\ ... np xy \ fa 2 ... p xy s 




K np xy ... na 2 y ) \p xy ... a 2 

where p xy = C ° v * y = is correlation coefficient of random quantities l x and l y (i.e. 

indicators of events from X). 

2.4 Polynomial distribution is a particular case of Binomial mul- 
tivariate distribution, when the latter is expressed by the 
partition of elementary events space 



£X x, 



When the generating N-set X consist of the events, which arising the partition Q = Y^, x 
then Binomial multivariate distribution of a random vector £ = (£ x ,x G X) is defined by 
the N parameters^ (n;p x ,x G X), where p x = P(x), YjxexPx = 1j an d it is Polynomial 
distribution with the given parameters. 

Hence, probabilities of the Binomial multivariate distribution, which is generated by the 
partition Q, is defined for every integer- valued vector h = (n x , x G X) from the simplex S N 
(because Yl x ex n % = n ) by the same formula as probabilities of corresponding Polynomial 
distribution 

bn(n;p x ,x G X) = P(| = n) = P(£ x = n x ,x G X) = — JJ[pJ 

11 



n x- xex 



x€X 



Since ^2 x pxPx = lj there are only N — 1 independent probabilities among N. 
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2.5 Binomial iV-variate distribution, which is generated by the 
set X, defines Polynomial 2 7V -variate distribution, which is 
generated by the set of terrace-events {ter(X), X C X}, but 
not visa versa 

Multivariate (N -variate) Bernoulli testing scheme of n tests with the generating set of 
events X, which obey eventological distribution {p(X), X C X}, defines N random quan- 
tities 

n 

each of them has Binomial distribution with parameters (n;p x = P(x)) and all together 
forms iV- variate random vector £ = (£ x , x G X), which is distributed by the Binomial 
multivariate (N -variate) law with 2 N parameters (n; {p(X), 0^IC X}), which contains 
amount of tests n and 2^ — 1 probabilities from eventological distribution of the generating 
set of events X (in other words, all 2^ probabilities p(X) without p(0))- 

The same Bernoulli multivariate testing scheme of n tests defines 2 N random quantities 

n 

ax)(uj) = X)W«)M' x ® ^ xW ' x e x, 

1=1 

each of them has Binomial distribution with parameters (n;p(X) = P(ter(X))), and all 
together forms 2^- variate random vector £ = (£ (X) , X C X) , which is distributed by the 
Polynomial multivariate (2 N -variate) law, which is generated by the set of terrace-events 
{ter(X), 

X C X} and which is defined by 2 N + 1 parameters (n; {p(X), X C X}. Those parameters 
contains amount of tests n and all 2 N probabilities p(X) from eventological distribution of 
the generating set X of events. 

Probabilities of the given Binomial and Polynomial multivariate distributions are bound 
for every ^-variate collections of nonnegative numbers n = {n x ,x G X} G [0,71]^ by the 
formula P(£ = n) = ^n-P(£ = ^) wnere summation is made by the all 2^- variate sets 
of nonnegative numbers n = (n(X), X C X) G S 2N for which n = J2xcx n ( X ) an< ^ a ^ so 
n x = J2 xe X n ( X )> ^ e X. 

Remark. For any Binomial iV-variate, which is generated by the set of events X, there 
is unique Polynomial 2 7V -variate distribution, which is generated by the set of terrace- 
events {ter(X),X C X}. The contrary is not true, i.e. for arbitrary Polynomial 2^- variate 
distribution, which is generated by the 2 N -set of events (those events form partition of Q), 
there are, generally speaking, (2^)! Binomial N- variate distributions, which is generated 
by the iV-sets of terrace-events X (appreciably depends on the way of partition events' 
labelling as subsets X C X and total amount of such ways is equal to (2^)!). 

3 Poisson multivariate distribution 

Poisson multivariate distribution is a discrete distribution of probabilities of a random 
vector £ = (£ X) x G X), which have values n = (n x , x G X) with the probabilities 

P(| = n) = P(& = n x , x G X) = 7r fl (Ap0, ^ X C X) = 
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where summation is made by collections of such nonnegative integer-valued numbers n = 
(n(X), 

7^ X C X) , for which there are iV equations = XLex n (^0> x £ and {A(X), 7^ 
X C X} with parameters: A(X) is an average number of coming of the terrace-event 



ter(X) = p| x p| 



x\ 



xex xex c 



in other words, average number of coming of the all events from X and there are not events 
from X c . 

A = J2 APO 

is an average number of coming at least one event from X, in other words, average number 
of coming of event x (union of all events from X). 

For example, when n — (0, . . . , 0) then 

P(|=(0,...,0))=P(^ = 0, xeX)=e~ x , 

when n — (0, . . . , 0, n x , 0, . . . , 0), x G X then 

P{£= (0,... ,0,^,0,. ..,0)) =P(£ x = n x , Z y = 0,y^x)=e- x [X ^ ]n{ * ] . 

If the vector n has one fixed component n x and other items are arbitrary: 

h = (-,...,-, n x , •,...,•), x G X, 

then 

p (l = (•,•••, •,...,•)) =P(C* = w x ) = e" 



is a formula of Poisson one-variate distribution with parameter X x of the random quantity 
£ x , where A^ = XLex ^(^) is defined for each x G X by the parameter of the Poisson 
one-variate distribution. 



3.1 Eventological interpretation 

Let there are countable sequence of n independent stochastic experiments. In the result of 
experiment n can ensue or not events from X. Possibilities p x = P(x) of events x G X are 
small, i.e. possibilities p(X), 7^ X C X of generated by them terrace-events ter(X), 7^ 
X C X are small too, and when n — > 00 then np(X) — > A(X), / I C I. Then random 
vector 

00 

I = (£,., x G X) = l x( „) , x G X 

n=l 

obeys multivariate (iV-variate) Poisson distribution with parameters {A(X), 0^IC X}. 
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Remark. It's incorrect to imagine that possibilities so tend to zero that only in n-th 
test np(X) = X(X), X C X. In truth, it's rather to believe that tending of possibilities to 
zero like that this equation is true for all first n tests. Thus, stochastic experiment consists 
in the sequence of n- series of independent tests (series of n tests), and this equation is 
true for all tests from n-series. Then n-series defines Binomial multivariate (iV-variate) 
distribution with parameters (n;p(X),0 ^ X C X), which by n — > 00 tends to the 
Poisson multivariate (iV-variate) with parameters (X(X), 0/ICl). 



3.2 Characteristics of the Poisson multivariate distribution 

Vector of mathematical expectations of the Poisson multivariate distribution is (E£ x , x G 
£) = (A x , x G X), where X x = XLex'H^)' x e Since Cov(^ 3 ,,^ J/ ) = X xy , where 

^xy = J2{ x ,y}CX 

{ x , y} Q 3£, so covariance matrix is equal to 

X X . . . X X y 

\^xy ' ' ' Xy 

In the case of two dimensions, when X = {x,y}, summation is making by the one 
parameter n(xy) = n({x,y}), which is changing in a Frechet bounds: 

P(| = n) = P(£ x = n x ,€ v = n y ) = 7v h (X(x), X(y), X(xy)) = 

= e _ x min ^ } [A(x)]"W [A(y)]"fa> [\(xy)] n W 
^ n(x)\ n(y)\ n(xy)\ 

where n(x) = n x — n(xy), n{y) = n y — n(xy), and A = X(x) + X(y) + X(xy). 

Vector of mathematical expectations Poisson two-variate distribution (E£ x , E£ y ) = 
(X x , X y ), where X x = X(x) + X{xy) and X y = X(y) + X(xy), covariance matrix is equal to 

/ X x X(xy)\ 
\X(xy) X y J ' 

because in the case of two dimensions X xy = X(xy). 



3.3 Poisson multivariate approximation 

If amount of independent experiments n is large-scale and possibilities p x = P(x) of events 
x G X is small (i.e. possibilities p{X), 7^ X C X of generated by them terrace-events 
ter(V),0 7^ X C X is small too) then for each collection of integer- valued numbers n = 
{n x , x G X) G [0,77,]-^ Binomial possibilities is expressed in the rough by terms of the 
Poisson multivariate distribution: 



b fl (n;p(X), 0^a), e -»E^**) J2 II [M fl\ (X) 
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where summation is applied to such sets (n(X),0 ^ X C X), for which n ^ 

and N equations n x = n(X), x E X are true. 

In the case of two dimensions, when X = {x,y}, summation is making by the one 
parameter n(xy) = n({x,y}), which is changing in so-called Frechet bounds: 

bn(n;p(x),p(y),p(xy)) « 

w e -n(p(*)+p(i,)+p(*i,)) min ^ nj/} [npfo)]"^) [np(y)} n ^ [np(xy)} n ^ 

^ n(x)\ n(y)\ n(xy)\ 

n(xy)=0 V 1 yyl V yJ 

where n(x) = n x — n(xy), n(y) = n y — n(xy). 

Poisson theorem (multivariate case). Let p x ^ 0, x E X when n — > oo, and 
np(X) — > A(X) for all nonempty subsets ^ I C I as that. Then for any collection of 
integer-valued numbers n = (n x , x E X) E [0, n] N (when n — > oo) 

b h (n; p(X) (A(X) , ^ X C X) , 

where 

f X C X) — e-SwW^ JT 

ft 

is Poisson multivariate possibility, and summation is applied to such sets n, for which 

n z = X^ex n ( x )i x EX. 

Proof. Because when n is large while n ^ J2xcx n (X) is true for any fixed n(X), X C 
X, then summation in formulas of Binomial 

b h (n;p(X),0^XCX)=J2 m n ("I {p W , ^ ^ X}) 

ft 

and Poisson 

7r ft (A(X), ^ X C X) = e-^W £ J] 



[\(X)] n W 



n(X)\ 

ft v 7 

possibilities is applied to the same sets h, for which n x = Ylxex n (^0i x E X. 
Now we show Poisson approximation for Polynomial possibilities 



mfl (n;{p(X),XCX})= n ] J] [p(X)]«W. (1) 

11 n(X)\ X c X 



xcx 



It should be pointed out, that for any fixed n(X), XCX and sufficiently large n there are 
follow equations: 

m { n(0)Mz),Mx),z^xcx}) {n; {p(X),X C X}) _ p(Z) (n(0) + l) 



m( n (0)+i,„(z)-i,{ n (i),z/ia}) (n; {p(X) ,ia}) n(Z)p(0) 
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where Z C X. By multiplying and dividing numerator and denominator by n and in 

consideration of n ( ) +1 ~ i anc l p(0) « 1, where ~ signify approximate equality with 
precision up to n -1 , we obtain 

p(Z)(n(0) + l)n _np(Z) n(0) + 1 1 _ rap(Z) 
n(Z)p(0) n n(Z) n p{0) n {Z) 

By the data of the theorem np(Z) —>■ X(Z), therefore 

m (n(0),»(z),{n(x),z^cx}) (n; {p(X),X C X}) A(Z) 



m(n(0)+l,n(Z)-l,{n(X),Z^CJ}) («! 0?P0 ,^^1}) 7l(Z) 

When n(X) =O,0/ICI, then 

m (n(0) , Ol ... 1 o)(n; (p(I),K X}) = |p( )]» = ^ " E ^) = ( 



(2) 



1-A 

n 



where A = Excj^W- After finding the logarithm of both parts of the equation and 
factorizing into Maclaurin series^) we have 



in[m( n ( ) A ... jO )(ri; {p(X),X C X})] = n ■ In (l - ^ J 



\ 2 

-A- — 

2n 



When n is large we conclude that 

m( n ( 0)j o,...,o) (n; {p(X), A" C X}) « e" A . (3) 
By the sequentially applying equation (j2J) to approximation (j2D we come to 

[A(X)]" (X) 



"i(„( ),{„(x),xa}) (n; {f>(A"), X C X}) « e A [ [ 



n(X)! ' 

i.e. Poisson approximation of the Polynomial possibility (pQ), from where the assertion of 
the theorem follows directly. 



4 Conclusion 

The multivariate generalizations of Binomial and Poisson distribution offered in the paper 
allow to consider any structure of dependences of the fixed set of events in sequence of 
independent multivariate tests and include as very special variant for the Polynomial dis- 
tribution long time known in probability theory. Offered multivariate discrete distributions 
pay increase in set of parameters for this unique opportunity. 



fc=i 



11 



References 

[1] Vorobyev O.Yu. Eventology. — Krasnoyarsk. — Siberian Federal University. — 2007. 
- 435p.; ISBN 978-5-7638-0741-7 



Oleg Yu. Vorobyev 

Siberian Federal University 

Institute of Mathematics 

Krasnoyarsk State Trade Economic Institute 

mail: Academgorodok, P.O.Box 8699 

Krasnoyarsk 660036 

Russia 

e-mail: vorob@akadem.ru 



Lavrentiy S. Golovkov 

Siberian Federal University 

Institute of Mathematics 

Krasnoyarsk State Trade Economic Institute 

mail: ul. Lidy Prushinskoi 2 

Krasnoyarsk, 660075 

Russia 

e-mail: lavrentiy.golovkov@gmail.com 



url: http://eventology-theory.ru 



12 



